O 



X 



Neutral stability height correction for ocean winds 



Ross N. Hoffman* 

O July 8, 2011 

(N 



Abstract 



I— s 

t> 

Adjusting ocean wind observations to a standard height, usually 10 m, requires 

the use of a boundary layer model, and knowledge of the thermodynamical variables. 
Height adjustment is complicated by the fact that a necessary parameter, the rough- 
ness height, cannot be given in a closed form solution. If only the wind and reporting 
height are known, the best that can be done is to assume neutral stability. The deter- 



mination of roughness height is analyzed and a simple approximation (used by Atlas 
et al.||2011 ) is derived in detail. This approximation is accurate for winds in the range 
of 1 — 30 m s 1 for neutral stratification and would be an excellent initial estimate for 
a Newton iteration to determine the roughness height precisely, whether or not neutral 
stability is assumed. 

> 

^2 1. Introduction 

Adjusting ocean wind observations to a standard height, usually 10 m, requires the use of a bound- 
ary layer model, and knowledge of the thermodynamical variables. Whichever PBL model is used, 
an exact solution requires iterating the constant flux layer equations. This is due to the fact that zo, 
^^ the roughness length, is an implicit function of the model variables over the oceans. The Charnock 

formula states that over the ocean zo and the surface stress magnitude |t| are linearly related by 



a , 

Z0= — T- (1) 

P8 

Here a = 0.0185 is the Charnock constant, g = 9.81 m s~ 2 , and p, the density of air, is assumed to 
be constant for the range of heights considered and equal to its surface value. Previously a = 0.032 
was the accepted value for the Charnock constant, but now it is known the Charnock constant is not 
actually a constant, but depends on the sea state, through the wave age ( Wu||1985 ). Accordingly, 



the Charnock constant is determined by the wave model in the ECMWF system, and normally is 
in the range 0.01 < a < 0.04, corresponding to sea states from swell to steep young ocean waves 
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( |Hersbach 2011). While values of a as large as 0.1 sometimes occur in the ECMWF system, a 
typical value is a = 0.018, which agrees with the value of 0.0185 used here and by Wu ( 1985). 
Note that Eq. Q neglects the contribution of molecular viscosity which is important at low wind 
speeds (e.g., Hersba ch|2011 ). However, at low wind speeds the height correction and consequently 
errors made in the height correction should be small. 
Surface stress determined from 

T = -pC d \V\V (2) 

also depends on zo through the neutral drag coefficient, and in the unstable case through the sim- 
ilarity function, usually denoted f(Ri) where Ri is the Richardson number. In Eq. ^ V is the 
vector wind at some height z, \V\ is the magnitude of the vector wind, the wind and stress vectors 
are assumed to be parallel for the range of heights considered, and the drag coefficient Cj is given 
by the product of the similarity function f(Ri) and the neutral drag coefficient 
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(3) 



Here the von Karman constant k = 0.4. See Hoffman and Louis ( 1990 1 for details 



NWP models usually "cheat" and use the value of x of the previous time step to find zo through 
the Charnock formula. Actually using old values to evaluate the dissipative terms can be a good 
policy as this can reduce computational instability. But outside of a model we must calculate zo 
implicitly. For this purpose we substitute the absolute value of Eq. ([2]) into Eq. ([T]) and then use the 
expresseion for Q/ to obtain 



zo 



-c d \v\ 
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k\V\ 
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= h(z ;\V\,z,Ri). 



(4) 



For neutral stratification, Ri — 0, and f(Ri) = 1. To solve Eq. (Eh we must iterate. To begin the 
process, Hoffman and Louis ( 1990) estimated Cd as a linear function of |V| and then obtained the 



initial estimate of zo from the Charnock relationship. Then Eq. ([4]), zo = h, is iterated. This con- 
verges to a good approximation within a few iterations. It is then possible to switch to a Newton 
iteration to solve ZQ — h = /(zo) = 0. The Newton method requires the partial derivative of h with 
respect to zo- This can be evaluated using the tangent linear code corresponding to the calculation 
of h by setting all inputs to zero except for that corresponding to zo, which is set to unity. The ad- 
vantage of the Newton method is that it iterates to machine precision in 2-4 steps from a reasonable 
start. (With an unreasonable start it can diverge.) With a solution exact to machine precision one 
can then skip the iteration in the adjoint and/or tangent models. As an alternative within the context 



of the ECMWF system, Hersbach (2011) describes an accurate fit for Q„ and zo as functions of 
neutral wind speed and the Charnock value, two parameters available in the interface between the 
ECMWF atmospheric and wave models. 



2. Height correction for ocean winds 

Knowing zq is equivalent to knowing the stress, and we can then solve Eq. ^ for |V| at any height, 



z, if only we know the Richardson number. In particular, Eq. (4 ) states that CJ |V| is conserved as 
we vary the height z. However Ri depends on knowing the stratification of the boundary layer. In 
what follows we assume that only the wind and reporting height are known. Then the best that can 
be done is to assume neutral stability, and in the rest of this treatment V will denote the neutral 

stability wind. Now, Cj \V\ is conserved, allowing us to determine the 10 m neutral wind speed 
| Vio| from an observation at some other height according to: 



[Viol 



log(lOAo) 

log (z/z ) 



IVI. 



(5) 



Note that according to Eq. ([5]), the ratio between neutral stability winds at two levels is entirely 
determined by zo and the two heights. 

Once zo is determined, the neutral wind, defined by 



T = -pC dn \V\V, 
is easily determined from knowledge of zo alone according to 



V 
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(6) 



(7) 



which is obtained by combining Eq. ([6]) with Eq. ([T]) and making use of Eq. ([3]). A few sample 
calculations using Eq. ([TJ) are presented in Table [T] for heights of 4, 10, and 19.5 m. In Table [T] 
we see that the variation in zo is two orders of magnitude greater than the variation in wind speed. 
Over this range of wind speed, the correction factors for determining |Vio| var y by as much as 5%. 
This variation is the same order of magnitude as the corrections, and is therefore worth accounting 
for. 



Table 1: Sample calculations based on Eq. §7j. The zo values are equal to 2~ J but have been 
multiplied by 10 6 for presentation in this table. The |Vio| values are in m s~ ! and are calculated 
using Eq. ([7j). The ratios in columns 4 and 5 are equal to the term in square brackets in Eq. ([5]). 



The last column contains the estimated value of zq from Eq. (12). 



j zoxlO 6 



14 
11 

8 



61 

488 

3906 



[Viol |Vi /V 4 | IV10/V19.5I 



5.4 
12.6 

28.2 



1.08 
1.10 
1.13 



0.947 
0.937 
0.922 



zo x 10 6 



68 

493 
3216 



3. Calculation of zo under neutral conditions 

To apply Eq. Q we still need to determine zq. Here we demonstrate a simple approximation. The 
motivation is that under neutral conditions, for some fixed height, we expect wind speed, surface 
stress, and roughness height to all increase together. Differentiating Eq. ([7]), we obtain 

d\V\ ( g \i _i 

zo 2 



dzo \ak 2 
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(8) 



Thus d\V\/dzo/ > provided z > e 2 zo- This holds for wind speeds less than hurricane strength 
and heights of several meters or more. The suggestion then is that zq should be a monotonically 
increasing function of |V|, and interpolation into a look-up table, or a simple fit to a set of exact 
values should work. 

For this investigation it is convenient to define 

y = log(z/z ) so that zo = ze~ y . (9) 

Then the square of Eq. ([7]) may be written as 

y 2 e -y = r= al l v 2 (10) 

gz 



Usually we will know V and z and hence y from Eq. (10). From y we then determine y and 



finally zo from Eq. ([9]). To determine y from y we tabulate or model y as a function of y, based 



on data obtained by calculating y from Eq. (10) for different values of y. Note that the values 



of the regression coefficients determined below are independent of the value of z or any of the 
other parameters, including the Charnock constant a. However, to create a relevant sample of y- 
values for fitting, we take z = 10 m, and vary zo- Below, as in Table [Tj we take values of zo evenly 
distributed in log space given by 2~ J for integer values of j. 

Figure [T^ plots y as a function of log(y) for j = 0, . . . , 30. Clearly a linear fit will work well 



over most the range of y. This is not unexpected since according to Eq. ( 10 ) log y = — y + 2 logy. 
As j increases, zo, |^io|» an d Y decrease, while y increases. For example, for j = 0, zo = 1 m and 
| Vio| = 134 ms _1 , while for j = 30, zo ~ 10~ 9 mand |Vio| = 0.04 m s _1 . Fitting points for j > 6, 
corresponding to |Vio| < 40 m s _1 (where the extreme point included is marked by the vertical line 
in the plots) we find that 

y = c + cilog(y) (11) 

with co = 3.7 and c\ = —1.165. 

Combining Eq. ([9]), Eq. ( fT0| ), and Eq. ( fTTj ), we obtain our estimate of zo 

zo = zexp[-(c + cilog[y])] =zexp[-(c + cilog[(afc 2 V 2 )/(sz)])]. (12) 

Figure [TJ? shows the error of the fit in log space. Values calculated using Eq. ( [12] ) are shown in the 
last column of Table [TJ for the cases listed. The differences are not tiny, but when we recalculate 
the ratios of the wind speeds in columns 4 and 5 of the table, the results are nearly the same. Using 
the same precision as in the table, the values are the same except that the value for IV10/V19 ,s\ is 
0.923. 



Figure 1: Fitting y = \og(z/zo) as a function of 7= (ak 2 / gz)V 2 . (a) The linear fit to values of y 
and log 7 for j — 6, . . . , 30 is plotted as a dotted line. The data values are plotted as dots, (b) Log 



residuals for Eq. ( 12 ) fit of zq. For this calculation we define the true values of zo and then y, 7 and 
£0 for the chosen value of z (10 m), using Eq. Q and Eq. ( fT2| ). The vertical lines identify j = 6. 
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4. Concluding remarks 



An approximation for roughness height (Eq. ( 12 )) is given that is accurate for winds in the range of 
1 — 30 m s _1 , assuming neutral stability and neglecting effects of molecular viscosity. (Also, note 
that coherent structures in the atmospheric boundary layer are not explicitly included in the simi- 
larity theory employed here.) Values of zo that vary with wind speed should be used in correcting 
ocean winds to a standard height. Typical corrections are in the range of 5-15%, so assuming a 
single value for zq will incur errors of a few per cent. Our approximation is also an excellent initial 
estimate to begin a Newton iteration to determine the roughness height precisely, whether or not 
neutral stability is assumed. 

In practice, the approximation derived here is adequate because errors due to other approxima- 
tions and assumptions are graver. For example, observations and meta-data associated with ship 
reports are often limited: information required to estimate atmospheric boundary layer stability 
may be lacking, anemometer heights may be unknown or incorrect, and effects due to ship motion 
and flow anomalies due to superstructure may not be accounted for. 
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